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We update the constraints on possible features in the primordial inflationary density perturbation 
spectrum by using the latest data from the WMAP7 and ACT Cosmic Microwave Background 
experiments. The inclusion of new data significantly improves the constraints with respect to older 
work, especially to smaller angular scales. While we found no clear statistical evidence in the data for 
extensions to the simplest, featureless, inflationary model, models with a step provide a significantly 
better fit than standard featureless power-law spectra. We show that the possibility of a step in the 
inflationary potential like the one preferred by current data will soon be tested by the forthcoming 
temperature and polarization data from the Planck satellite mission. 

PACS numbers: 98.80.Cq, 98.70.Vc, 98.80.Es 
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I. INTRODUCTION 



Current cosniological observations can be explained 
in terms of the so-called concordance ACDM model in 
which the primordial fluctuations are created during an 
early period of inflationary expansion of the Universe. 
In particular, the spectrum of anisotropics of the cosmic 
microvifave background (CMB) is in excellent agreement 
with the inflationary prediction of adiabatic primordial 
perturbations with a nearly scale-invariant power spec- 
trum [lH5|. In its simplest implementation, inflation is 
driven by the potential energy of a single scalar field, the 
inflaton, slowly rolling down towards a minimum of its 
potential; curvature perturbations, that constitute the 
primordial seeds for structure formation, are originated 
during the slow roll from quantum fluctuations in the 
inflaton itself. The scale invariance of the spectrum is 
directly related to the flatness and smoothness of the in- 
flaton potential, that are necessary to ensure that the 
slow-roll phase lasts long enough to solve the paradoxes 
of the Big Bang model. 

However, in more general inflationary models, there 
is the possibility that slow roll is briefly violated. This 
naturally happens in theories with many interacting 
scalar fields, as it is the case, for example, in a class 
of multifield, supergravity-inspired models [a, 0|) where 
supersymmetry-breaking phase transitions occur during 
inflation. These phase transitions correspond to sudden 
changes in the inflaton effective mass and can be mod- 
eled as steps in the inflationary potential. If the transi- 
tion is very strong, it can stop the inflationary phase as 
it happens in the usual hybrid inflation scenario; on the 
contrary, inflation can continue but the inflationary per- 
turbations and thus the shape of the primordial power 
spectrum are affected. Departures from the standard 
power-law behaviour can also be caused by changes in 
the initial conditions due to trans-planckian physics ^- 
Il0 | or to unusual initial field dynamics [ill, II4I 



A violation of slow-roll will possibly lead to detectable 
effects on the cosmological observables, or at least to the 
opportunity to constraint these models by the absence 
of such effects. In particular, step-like features in the 
primordial power spectrum have been shown |13l . Il4| to 
lead to characteristic localized oscillations in the power 
spectrum of the primordial curvature perturbation. Such 
oscillations have been considered as a possible explana- 
tion to the "glitches" observed by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) in the temperature 
anisotropy spectrum of the CMB, although the WMAP 
team notes that these could be caused simply by having 
neglected beam asymmetry, the gravitational lensing of 
the CMB, non-gaussianity in the CMB maps and other 
"small" (< 1%) contributions to the covariance matrix. 
In the following we will assume that these features have 
indeed a cosmological origin as in the class of extended 
models described above and we will use CMB data to con- 
strain the phenomenological parameters describing the 
step in the inflaton potential. 

Constraints on oscillation in the primordial perturba- 
tion spectrum, as well as best-fit values for the step pa- 
rameters, have been previously derived in Refs. [15i - tlil |. 
Here we improve on the previous analyses in several as- 
pects. First, we use more recent CMB data, in particular 
the WMAP 7-year and the Atacama Cosmology Tele- 
scope (ACT) data. This allows us to derive tighter con- 
straints on the parameters; in particular we get an upper 
limit on the step height (related to the amplitude of os- 
cillations) that is independent on the position of the step 
itself in the prior range considered. We also find a clear 
correlation between the position and the height of the 
step. Secondly, we generate mock data corresponding to 
the model providing the best-fit to the WMAP data, and 
use these data to assess the ability of the Planck satellite 
to detect the presence of oscillations in the primordial 
spectrum. 

The paper is organized as follows: in Section II we de- 
scribe the evolution of perturbations in interrupted slow 



roll and the phenomenological model used to describe 
a step in the inflationary potential. In Section III we 
discuss the analysis method adopted. In Section IV we 
present the results and in Section V we derive our con- 
clusions. 



II. INFLATIONARY PERTURBATIONS IN 
MODELS WITH INTERRUPTED SLOW ROLL 

Steps in the potential can naturally appear in "multi- 
ple inflation" models, where the inflaton field (j) is gravi- 
tationally coupled to a "flat direction" field p (belonging 
to the visible sector of the theory), i.e. a direction in field 
space along which the potential vanishes. The p-field can 
undergo a symmetry-breaking phase transition and ac- 
quire a vacuum expectation value (p). The gravitational 
coupling between the p and the infiaton field will cause 
the effective mass-squared of the latter to change; for ex- 
ample, in the case in which the coupling between the two 
fields is described by a term A(/)^p^/2 in the Lagrangian, 
the inflaton mass-squared after the phase transition will 
become m^g((/)) = ttiq -I- A(p^). It is worth noticing that 
the presence of fiat field directions also opens the possi- 
bility to have inflation with a curved trajectory in field 
space; however, in the following, we will disregard this 
scenario. 

The exact behaviour of the infiaton mass will depend 
by the dynamics of the phase transition; however, this 
is so fast that the p-field reaches the minimum of its 
potential very rapidly. It is then very reasonable to model 
the infiaton mass in a phenomenological way as 
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1 + c tanh 



Here, the parameter b is of the order of the critical value 
of the infiaton field for which the phase transition occurs, 
c is the height of the step (related to the change in the 
inflaton mass) and d is its width (related to the duration 
of the phase transition). In the following we shall work 
in reduced Planck units {c = h — SttG = 1), so that all 
dimensional quantities like rn, b and d should be multi- 
plied by the reduced Planck mass Mp = 2.435 x 10^^ GeV 
in order to get their values in physical units. 

Let us now briefly recall how to compute the spectrum 
of primordial perturbations, as discussed in details by 
Adams et al. [ij]. For the moment, we do not specify 
the exact form of the inflaton potential V{4)); we will 
return on this in the next Section. In the case of scalar 
perturbations, it is useful to define the gauge-invariant 
quantity |2C| u = ~zTZ, where z — a(j>/H, a is the scale 
factor, H is the Hubble parameter, TZ is the curvature 
perturbation, and dots denote derivatives with respect 
to the cosmological time t. The Fourier components of u 
evolve according to: 



^k^ \k ) Ufc = , 



where k is the wavenumber of the mode, and primes de- 
note derivatives with respect to conformal time 77. When 
k'^ ^ z" /z, the solution to the above equation tends to 
the free-field solution Uk = e~*'"'''/-\/2fc. 

In the slow-roll approximation, z" /z ~ 2a^H^. How- 
ever, in the models considered here this expectation can 
be grossly violated near the phase transition, and the 
time evolution of z has to be derived by solving the equa- 
tions for the inflaton field and for the Hubble parameter: 
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F(0). 
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Once the form of the potential is given, these can be 
integrated to get H and (/>, and thus z, as a function of 
time. At this point, it is possible to integrate Eq. ([2]) 
to get Ukir]) for free-field initial conditions when /c^ ^ 
z" / z. Finally, knowing the solution for the mode fc, the 
power spectrum of the curvature perturbation Pji can be 
computed by means of 
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evaluated when the mode crosses the horizon. The re- 
sulting spectrum for models with a step in the potential 
is essentially a power-law with superimposed oscillations; 
thus, asymptotically, the spectrum will recover the famil- 
iar fc""'^^ form typical of slow-roll inflationary models. 

In practice, however, one has to relate the horizon size 
at the step with a physical wavenumber. For a general 
(]\ wavenumber ki, one can write fc^ = ai,Hi, — aonde~^*-ff*, 
where a^, and il* are the scale factor and the Hubble 
parameter at the time the mode crossed the horizon, 
flond is the scale factor at the end of inflation, and iV^ 
is the number of e-fold taking place after the mode left 
the horizon. We choose Ni, = 50 for the pivot wavenum- 
ber k^ = ko — 0.0025 Mpc~ . A different choice would 
correspond to a translation in the position of the step in 
(j) and would thus be highly degenerate with b. For this 
reason we do not treat iV^ as a free parameter, consistent 
with what has been done in previous studies ^16^ .17:]. 



(2) 



III. ANALYSIS METHOD 

We compare the theoretical predictions of a class of 
inflationary models with a step in the infiaton potential 
with observational data. We use a modified version of 
the CAMB code that solves Eqs. ©-(HI) numerically us- 
ing a Bulirsch-Stoer algorithm in order to compute the 
initial perturbation spectrum ([5]) and, from that, the 
CMB anisotropy spectrum for given values of the rele- 
vant parameters describing the model. CAMB is then in- 
terfaced with a modified version of the Markov Chain 
Monte Carlo (MCMC) package CosmoMC '21], that we 
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FIG. 1. Inflationary potential Q for m = 7.5 x lO"*^. The 
solid thin black line corresponds to a smooth (c = 0) chaotic 
potential m^cj)^ /2. The long-dashed red curve has b = 14.23, 
c = 0.001 and d — 0.025 and roughly corresponds to the 
spectrum giving the best fit to the WMAP7 data (see Sec. IIVI 
below). The other curves correspond to b — 14.23, c = 0.005, 
d = 0.025 (blue short-dashed), b = 14.23, c = 0.005, d = 0.05 
(green dotted) and b = 14.25, c = 0.005, d = 0.025 (magenta 
dot-dashed) . 



use to find the best-fit value of the parameters, to recon- 
struct their posterior probability density function, and to 
infer constraints on the parameter themselves. 

Models. We consider a chaotic inflation potential of 
the form F(0) = m^jj ((/))0^/2. Using Eq. ^, this corre- 
sponds to a potential 



V{cl>) = im^c/)^ 
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In Fig. [T] we show the shape of this potential for 
TO = 7.5 X 10~^ and different values of the step parame- 
ters (close to the best-fit values), compared to a smooth 
TO^(?!)^/2 potential (c = 0). 

The potential ^ uniquely defines the spectrum of per- 
turbations Vn- The parameters that define the primor- 
dial spectrum and the initial conditions for the evolution 
of cosmological perturbations are then the inflaton mass 
TO and the step parameters b, c and d. The infiaton mass 
sets the overall scale for the potential and consequently 
for the amplitude of the perturbations; it can then be 
traded, in the Monte Carlo analysis, for the more famil- 
iar parameter A^, i.e., the amplitude of the primordial 
spectrum at the pivot wavenumber feg = 0.0025 Mpc~ . 
On the other hand, as already noted above, a step in the 
potential produces a perturbation spectrum with oscil- 
lations superimposed over a smooth power law. In the 
case of the potential (j6l), the underlying power-law has 
a fixed spectral index Ug = 0.96. In Fig. [2] we show 
the primordal spectrum for different values of the step 
parameters. 

The results obtained in the case of a specific poten- 
tial will be, by definition, model-dependent. However, as 
argued in Ref. [17| , the issue of model dependence can 



FIG. 2. Primordial power spectrum for an inflationary poten- 
tial of the form ([6]). The values of the step parameters are the 
same as in Fig. ^, namely: 6 = 14.23, c = 0.001, d = 0.025 
(red long-dashed), 6 = 14.23, c = 0.005, d = 0.025 (blue 
short-dashed), b = 14.23, c = 0.005, d = 0.05 (green dotted) 
and b = 14.25, c = 0.005, d = 0.025 (magenta dot-dashed). 



be alleviated in a phenomenological way by restoring the 
spectral index as a free parameter, i.e., by defining the 

gen 



"generalized" spectrum Vj^ as 



n^"(fc)=n"(fc) 



-0.96 



(7) 



where 7^^(fc) is the spectrum induced by the chaotic po- 
tential ^. Since the latter has a overall tilt of 0.96, Us 
will describe the overall tilt of the generalized spectrum. 

Summarising, we consider two classes of models. Mod- 
els belonging to the first class (referred to as class A) cor- 
responding to the potential ([5]), are described by eight 
parameters: the physical baryon and cold dark matter 
densities ujb = Vlbh? and loc — ^ch?, the ratio 9 be- 
tween the sound horizon and the angular diameter dis- 
tance at decoupling, the optical depth to reionization r, 
the parameters 5, c and d of the step-infiation model, 
and the overall normalization of the primordial power 
spectrum As (equivalent to specifying m? as discussed 
above). Models in the second class, referred to as class 
B, correspond to the generalized step model ([7]) and are 
described by the effective tilt Ug in addition to the eight 
parameters of the first class. In both cases, we consider 
purely adiabatic initial conditions, impose flatness and 
neglect neutrino masses. We limit our analysis to scalar 
perturbations. 

Priors. Apart from the hard-coded priors of CosmoMC 
on Ho (40 km s^^ Mpc^^ < Hq <100 km s^^ Mpc^^) 
and the age of the Universe (10 Gyr < tg < 20 Gyr), 
we impose flat priors on Wh, Uc, 9, t and, when con- 
sidered. Us and a logarithmic prior on As- As we shall 
see, for these parameters the width of the posterior is 
much smaller than the prior range, so that the latter is 
not really relevant. For the step parameters, the situa- 
tion is complicated by the fact that the likelihood (and 



the posterior) does not go to zero in certain directions of 
the subspace. This happens in particular for very small 
values of c, for which the spectrum becomes indistin- 
guishable from a power law, and for values of b either too 
large or too small so that the features in the spectrum 
are moved outside the range of observable scales. Then 
we choose for b a flat prior 13 < 6 < 15, that roughly en- 
compasses said range. In the case of c and d, since we do 
not have any a priori information on these parameters, 
not even on their order of magnitude, we find convenient 
to consider a logarithmic prior on both of them. Hence, 
we take (in the following, log a; denotes the base 10 log- 
arithm) —6 < logc < — 1 and —2.5 < logd < —0.5. 
Additionally, since the combination c/d? is better con- 
strained by the data than d alone, we also impose a prior 
— 5 < log(c/d^) < 3. Finally, we recall that, since the 
posteriors for b, log c and log d do not necessarily vanish 
at the edge of the prior range, all integrals of the prob- 
ability density function depend on the extremes of inte- 
gration and are thus somewhat ill-defined. Care should 
then be taken when quoting confidence limits in the b, 
logc and logd subspaces. 

Datasets. We perform the statistical analysis for 
each of the models by comparing the theoretical predic- 
tions to two different datasets. The first includes the 
WMAP 7-year temperature and polarization anisotropy 
data (WMAP7). The likelihood is computed using the 
the WMAP likelihood code publicly available at the 
LAMBDA websiteQ. We marginalize over the amplitude 
of the Sunyaev-Zel'dovich signal. The second dataset in- 
cludes the WMAP7 data with the addition of the small- 
scale CMB temperature anisotropy data from the ACT 
experiment. For the ACT dataset we also consider two 
extra parameters accounting for the Poisson and clus- 
tering point sources foregrounds components. The ACT 
dataset is considered up to ^max = 2500. 

Other than deriving the limits on the models from 
existing data, we also assess the ability of future ex- 
periments, in particular of the Planck satellite, to im- 
prove these constraints. In order to do this, we simulate 
"mock" data corresponding to the step model that yields 
the best-fit to the WMAP 7 and then perform a statis- 
tical analysis on these data as if they were real. The 
forecast method we use is identical to the one presented 
in [22] and we refer to this paper for further details and 
references. The synthetic dataset is generated by consid- 
ering for each Q a noise spectrum given by: 



Experiment Channel FWHM AT/T AP/T 



Planck 


70 


14' 


4.7 


6.7 


/sky = 0.85 


100 


10' 


2.5 


4.0 




143 


7.1' 


2.2 


4.2 



TABLE I. Planck [231 experimental specifications. Chan- 
nel frequency is given in GHz, FWHM in arcminutes and 
noise per pixel for the Stokes I (AT/T), Q and U parameters 
(AP/T) is in [10VK/K], where T = Tcmb = 2.725K. In the 
analysis, we assume that beam uncertainties and foreground 
uncertainties are smaller than the statistical errors. 




FIG. 3. Model likelihood as a function of b for model A (thin 
curves) and B (thick curves) using WMAP7 data (dashed 
curves) and the WMAP7-I-ACT dataset (solid curves) 



sensitivity a hy w ^ = (^bcr)^. 
eters are reported in Table |I| 



The experimental param- 



Together with the primary anisotropy signal we also 
take into account information from CMB weak lensing, 
considering the power spectrum of the deflection field C^'^ 
and its cross correlation with temperature maps Cf"^. 

Analysis. We derive our constraints from paral- 
lel chains generated using the Metropolis-Hastings algo- 
rithm. We use the Gelman and Rubin R parameter to 
evaluate the convergence of the chains, demanding that 
i? — 1 < 0.03. The one- and two-dimensional posteriors 
are derived by marginalizing over the other parameters. 



Ne ^ w-'^ cxp{£{e + l)8ln2/0l), 



(8) 



where 6'b is the full width at half maximum (FWHM) 
of the beam assuming a Gaussian profile and where w^^ 
is the experimental power noise related to the detectors 



^ http://lanibda.gsfc.nasa.gov/ 



TABLE II. Best-fit values for the parameters of the primordial 
spectrum. 





Model A 


Model A 


Model B 


Model B 


Parameter 


WMAP7 


WMAP7+ACT 


WMAP7 


W7+ACT 


b 


14.23 


14.25 


14.24 


14.25 


logc 


-3.11 


-2.71 


-2.97 


-2.67 


logd 


-1.58 


-1.60 


-1.65 


-1.45 


Us 


- 


- 


0.953 


0.959 


ln[10i°yl.] 


3.08 


3.06 


3.07 


3.08 


x' 


7469.4 


7489.6 


7467.9 


7491.4 




Model A, WMAP7 

Model B, WMAP7 

Model A, WMAP7+ACT 

Model B, WMAP7+ACT 




Model A, WMAP7 

Model B, WMAP7 

Model A, WMAP7+ACT 

Model B, WMAP7+ACT 



Model A, WMAP7 - 
Model B, WMAP7 ■ 
Model A, WMAP7+ACT 
Model B, WMAP7+ACT ■ 



Log[cl 



Log[d] 



FIG. 4. One-dimensional posterior probability density for b (left panel), logc (middle panel) and logd (right panel) for model 
A (thin curves) and B (thick curves) using WMAP7 data (dashed curves) and the WMAP7+ACT dataset (solid curves) 



IV. RESULTS AND DISCUSSION 

We first consider the WMAP7 and WMAP7-I-ACT 
datasets. We find that the ACDM fit to both datasets 
can be improved by the inclusion of a step in the infla- 
tionary potential, in both cases when the scalar spectral 
index is being fixed to n^ — 0.96 (Model A), and when it 
is being treated as a free parameter (Model B). The best- 
fit values for the step parameters are reported in Table 
im We also show the full likelihood for b in Fig. |31 It 
can be seen that in all cases the maximum in the likeli- 
hood occurs for b ~ 14.2; as we show below, this is due 
to oscillations placed in correspondence to the WMAP 
glitches at £ ~ 20 and £ ~ 40 and thus able to improve, 
for suitable values of the other parameters, the goodness- 
of-fi-t with respect to the vanilla ACDM model. We found 
that in case of the WMAP7 analysis the best-fit vanilla 
A-CDM model is at about Axlg ~ 6 from the global best 
fit with features. 

As long as bayesian statistics is concerned, the actual 
probability density distribution for a parameter is not 
given by the likelihood (the probability of the data given 
the parameters) but instead by the posterior (the prob- 
ability of the parameters given the data). In Fig. |4] we 
show the one-dimensional posterior distributions for the 
step parameters b, logc and logd. It can be noted that 
the posterior for b has a peculiar shape, presenting a peak 
for b ~ 14.2 and a fairly wide dip for 6 ^ 14. The peak 
traces the peak in the likelihood discussed above. The 
decrease for 5 < 14 is instead due to the fact that, low- 
ering 6, the oscillations are moved to larger multipoles 
where they tend to spoil the ACDM fit unless c is set to 
a very small value. 

This is clearly illustrated in Fig. [SJ where we compare 
the WMAP7 data with three realizations of the CMB 
spectrum: the ACDM best fit to the WMAP data, the 
generalized step model best fit to the same data (corre- 
sponding to the third column of Table |ll]) , and a gener- 
alized step model with the same parameters as the best 
fit, with the exception of b that is set to 5 = 13.9. It is 
clear, especially from the second panel, that for b — 14.2 
the oscillations improve the fit in the region 20 < ^ < 50. 
On the other hand, when b = 13.9 the height of the first 
peak is diminished so that the predicted spectrum is com- 



pletely at variance with the data. The posterior does not 
drop to zero because it still exist a fair amount of param- 
eter space, i.e., models with low c, than can fit the data 
even with the oscillations placed in the "wrong" place. 
The posterior going to a constant value at the edges of 
the prior range is instead related to the oscillations being 
moved out of the observable scales. The inclusion of the 
ACT data in addition to WMAP7 helps in constraining 
small values of 5, i.e., oscillations at small scales (large 
ts). 

The shape of the log c posterior is typical of a quantity 
parametrizing the amplitude of a non-standard effect: it 
is constant for "small" values of the parameter (when 
the step model becomes indistinguishable from standard 
ACDM), and then rapidly vanishes above a critical value. 
It can be seen that the probability density becomes half 
of its asymptotic value at c = for c > 10^^. Finally, 
the posterior for logc clearly shows that this parameter 
is largely unconstrained by data. 

We do not quote one-dimensional confidence limits on 
the parameters because, as noted in Sec. IIIIl the pos- 
teriors do not vanish at the edge of the prior range and 
in this case the confidence limits depend on the integra- 
tion range chosen. However, for illustrative purposes, in 
Fig. [B]we show the 2-dimensional 95% confidence regions, 
computed assuming that the posterior vanishes outside 
the prior range, in the (b — log c) plane. It is clear from 
the plots that there is a region below 6 = 14 where the 
data are more sensitive to the value of c; this is related as 
noted above to the oscillations being placed in the region 
where the data are more accurate and favour a smooth 
spectrum over one with oscillations. 

The results presented here are fully compatible with 
the analysis made by [1^ where the WMAP 5 dataset was 
considered. The apparently different value for the best- 
fit b parameter found in that paper is due to the differ- 
ent choice of the pivot scale (/cq = 0.05 Mpc^^ instead of 
ko = 0.0025 Mpc~^ as assumed in our analysis). We have 
checked that performing the analysis on the WMAP7 
dataset with the assumption of ko = 0.05 Mpc~^ results 
in a best-fit value of 5 ^ 14.7 in agreement with the re- 
suhs of [3. 

Finally, we show our results on the sensitivity of Planck 
to the step parameters. We have assumed as a fidu- 
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FIG. 5. Upper panel: CMB anisotropy spectrum for the 
ACDM (red solid line) and generalized step model (blue long 
dashed line) best fits, and for a step model with 6 — 13.9 
(black short dashed line), compared with the WMAP7 data. 
Lower panel: Zoom of the region £ < 60, showing the im- 
proved fit of the step model. 



TABLE in. Parameter constraints from Planck. 



Parameter 



Model A 



Model B 



b 


14.200 ±0.010 


14.200 ±0.011 


logc 


-3.00 ±0.32 


-3.00 ±0.34 


logd 


-1.66 ±0.22 


-1.64 ±0.23 


ris 


0.96 (fixed) 


0.957 ±0.007 


10" AJ 


3.073 ± 0.016 


3.074 ± 0.016 



cial model a generalized step model with b — 14.2, 
logc = -2.97, logd = 1.65, Us = 0.953, As = 2.16 x lO'^ 
(basically corresponding to the Model B best-fit to the 
WMAP7 data, i.e, the third column of Tab. El)- The one- 
dimensional posteriors for b, logc and logd are shown in 
Fig. [71 while in Tab. IIIII we report the mean values for 
the primordial spectrum parameters together with their 
2(7 error. As we can see the prior range dependence goes 
away with Planck data and we can quote marginalised 
credible intervals. We also show the two-dimensional pos- 
teriors for the step parameters in Fig. [8l It is evident 
that the Planck data will greatly increase the precision 
to which the step parameters can be measured; in par- 
ticular, a detection of oscillations will be possible. 
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FIG. 6. 95% two-dimensional confidence region in the (6- 
logc) plane. The four panels correspond to, from left to right 
and from top to bottom: class A, WMAP7±ACT; class B, 
WMAP7+ACT; class A, WMAP7; class B, WMAP7. 



CONCLUDING REMARKS 



We have considered inflation models with a small- 
amplitude step-like feature in the inflaton potential. Fea- 
tures of these kind can be due for example to phase 
transitions occurring during the slow roll in multi-field 
inflationary models. In these models the primordial per- 
turbation spectrum has the form of a power-law (as in 
the standard featureless case) with superimposed oscilla- 
tions, localized in a finite range of scales that basically 
depends on the position of the step in the potential. We 
have compared the theoretical predictions of a specific 
model, i.e., chaotic inflation, and of a more general phc- 
nomenological model to the WMAP7 and ACT data, in 
order to find constraints on the parameter describing the 
model. We have also studied the possibility of detecting 
the oscillations with the upcoming Planck data in the 
case they really exist. 

We have found that models with features can improve 
the fit to the WMAP7 data when the step in the potential 
is placed in such a way as to produce oscillations in the 
region 20 ^ £ < 60, where the WMAP7 data shows some 
glitches. We found no further evidence for small scales 
glitches from the recent ACT data, this is fully consistent 
with the recent analysis of Q. We have also found that 
models with too high a step are excluded by the data. 
Finally, assuming as a fiducial model the generalized step 
model that provides the best fit to the WMAP7 data, we 
have found that the Planck data will allow to measure 
the parameters of the model with remarkable precision, 
possibly confirming the presence of glitches in the region 
20 < ^ < 60. 






FIG. 7. One-dimensional posterior probability density for b (left panel), logc (middle panel) and logd (right panel) derived 
from the mock Planck data, for models of class A (dashed curves) and B (solid curves). 
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